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We present the first observation and cross section measurement of exciusive dijet production in 
pp interactions, pp — > p + dijet + p. Using a data sampie of 310 pb" 1 collected by the Run II Collider 
Detector at Fermilab at y / s=1.96 TeV, exclusive cross sections for events with two jets of transverse 
energy E^ 1 > 10 GeV have been measured as a function of minimum E^* ■ The exclusive signal is 
extracted from fits to data distributions based on Monte Carlo simulations of expected dijet signal 
and background shapes. The simulated background distribution shapes are checked in a study of a 
largely independent data sample of 200 pb _1 of 6-tagged jet events, where exclusive dijet production 
is expected to be suppressed by the J z = total angular momentum selection rule. Results obtained 
are compared with theoretical expectations, and implications for exclusive Higgs boson production 
at the pp Large Hadron Collider at y^s =14 TeV are discussed. 
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I. INTRODUCTION 

Exclusive dijet production in pp collisions is a process 
in which both the antiproton and proton escape the in- 
teraction point intact and a two-jet system is centrally 
produced: 

P + P^P' + {jet 1 +jet 2 )+p'. (1) 

This process is a particular case of dijet production in 
double Pomeron exchange (DPE), a diffractive process in 
which the antiproton and proton suffer a small fractional 
momentum loss, and a system X containing the jets is 
produced, 

P + P - [p 1 + Pp] + [p' + P p ] ^p' + X+p', (2) 

where P designates a Pomeron, defined as an exchange 
consisting of a colorless combination of gluons and/or 
quarks carrying the quantum numbers of the vacuum. 

In a particle-like Pomeron picture (e.g. see [l|), the 
system X may be thought of as being produced by the 
collision of two Pomerons, Pp and P p , 

Pp + P p ^ X Y P/ p + (jeh + jeh) + Y P/p , (3) 

where in addition to the jets the final state generally con- 
tains Pomeron remnants designated by Yp/p and Yp/ P - 
Dijet production in DPE is a sub-process to dijet produc- 
tion in single diffraction (SD) dissociation, where only the 
antiproton (proton) survives while the proton (antipro- 
ton) dissociates. Schematic diagrams for SD and DPE 
dijet production are shown in Fig. [1] along with event 
topologies in pseudorapidity space (from Ref. [IJ). In 
SD, the escaping p is adjacent to a rapidity gap, defined 
as a region of pseudorapidity devoid of particles 0]. A 
rapidity gap arises because the Pomeron exchanged in a 
diffractive process is a colorless object of effective spin 
J > 1 and carries the quantum numbers of the vacuum. 
In DPE, two such rapidity gaps are present. 

Dijet production in DPE may occur as an exclusive 
process (J) with only the jets in the final state and no 
Pomeron remnants, either due to a fluctuation of the 
Pomeron remnants down to zero or with a much higher 
cross section in models in which the Pomeron is treated 
as a parton and the dijet system is produced in a 2 — ^ 2 
process analogous to 77 — > jet + jet @. 

In a special case exclusive dijets may be produced 
through _an intermediate state of a Higgs boson decay- 
ing into bb: 

Pp + P p ^H°^(b^ j eh) + {b^ j et 2 ). (4) 
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FIG. 1: Illustration of event topologies in pseudorapidity, 
rj, and associated Pomeron exchange diagrams for dijet pro- 
duction in (a) single diffraction and (b) double Pomeron ex- 
change. The shaded areas on the left side represent "underly- 
ing event" particles not associated with the jets [from Ref. 




FIG. 2: Leading order diagrams for (a) exclusive dijet and 
(b) exclusive Higgs boson production in pp collisions. 



Exclusive production may also occur through a t- 
channcl color-singlet two gluon exchange at leading order 
(LO) in perturbative quantum chromo-dynamics (QCD), 
as shown schematically in Fig. [5] (a), where one of the two 
gluons takes part in the hard scattering that produces the 
jets, while the other neutralizes the color flow [6j. A simi- 
lar diagram, Fig.[2](b), is used in Q to calculate exclusive 
Higgs boson production. 

Exclusive dijet production has never previously been 
observed in hadronic collisions. In addition to providing 
information on QCD aspects of vacuum quantum num- 
ber exchange, there is currently intense interest in using 
measured exclusive dijet production cross sections to cal- 
ibrate theoretical predictions for exclusive Higgs boson 
production at the Large Hadron Collider (LHC). Such 
predictions are generally hampered by large uncertain- 
ties due to non-perturbative suppression effects associ- 
ated with the rapidity gap survival probability. As these 
effects are common to exclusive dijet and Higgs boson 
production mechanisms, dijet production potentially pro- 
vides a "standard candle" process against which to cali- 
brate the theoretical models [1, 0] ■ 

In Run I (1992-96) of the Fermilab Tevatron pp col- 
lider operating at 1.8 TeV, the Collider Detector at Fer- 
milab (CDF) collaboration made the first observation of 
dijet production by DPE) [2| using an inclusive sample 
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of SD events, pp — ► p'X, collected by triggering on a p 
detected in a forward Roman Pot Spectrometer (RPS). 
DPE dijet events were selected from this sample by re- 
quiring, in addition to the p detected by the RPS, the 
presence of two jets with transverse energy Et > 7 GeV 
and a rapidity gap in the outgoing proton direction in 
the range 2.4 < rj < 5.9 [8]. In the resulting sample of 
132 inclusive DPE dijet events, no evidence for exclusive 
dijet production was found, setting a 95 % confidence 
level upper limit of 3.7 nb on the exclusive production 
cross section. At that time, theoretical estimates of this 
cross section ranged from ~ 10 3 larger *9] to a few times 
smaller |6j than our measured upper bound. More data 
were clearly needed to observe an exclusive dijet signal 
and test theoretical predictions of kinematical properties 
and production rates. 

In Run II-A (2001-06), with the Tevatron providing 
pp collisions at y/s = 1.96 TeV, two high statistics data 
samples of DPE dijet events were collected by the up- 
graded CDF II detector: one of inclusive dijets, and an- 
other largely independent sample of &-quark jets. The 
analysis of these data is the subject of this paper. The 
results obtained provide the first evidence for exclusive 
dijet production in pp collisions. 

The paper is organized as follows. In Sec. [Ill we present 
the strategy employed to control the experimental is- 
sues involved in searching for an exclusive dijet signal. 
We then describe the detector (Sec. IIII[) . the data sam- 
ples and event selection fSec. IIV|) . the data analysis for 
inclusive DPE (Sec. [V} and exclusive dijet production 
(Sec. IVI[) . results and comparisons with theoretical pre- 
dictions (Sec. IVII[) . and background shape studies us- 
ing heavy flavor quark jets (Sec. IVIIip . Implications for 
exclusive Higgs boson production at the Large Hadron 
Collider are discussed in Sec. IIX1 and conclusions are 
presented in Sec. [Xj 



II. STRATEGY 

Exclusive dijet production is characterized by two jets 
in the final state and no additional final state particles 
except for the escaping forward proton and antiproton. 
Therefore, searching for exclusive dijet production would 
ideally require a full acceptance detector in which all fi- 
nal state particles are detected, their vector momenta are 
measured, the correct particles are assigned to each jet, 
and "exclusivity" is certified by the absence of any addi- 
tional final state particle(s). Assigning particles to a jet 
is a formidable challenge because the detector threshold 
settings used to reduce noise may inadvertently either 
eliminate particles with energies below threshold or else 
result in noise being counted as additional particles if the 
thresholds are set too low. To meet this challenge, we de- 
veloped a strategy incorporating detector design, online 
triggers, data sets used for background estimates, and 
an analysis technique sensitive to an exclusive signal but 
relatively immune to the above mentioned effects. 



The exclusive signal is extracted using the "dijet mass 
fraction" method developed in our Run I data analysis. 
From the energies and momenta of the jets in an event, 
the ratio Rjj = Mjj/Mx of the dijet mass Mjj to the 
total mass Mx of the final state (excluding the p and 
p) is formed and used to discriminate between the sig- 
nal of exclusive dijets, expected to appear at Rjj = 1, 
and the background of inclusive DPE dijets, expected 
to have a continuous distribution concentrated at lower 
Rjj values. Because of smearing effects in the measure- 
ment of Ej? and rf et and gluon radiation from the jets 
the exclusive dijet peak is broadened and shifts to lower 
Rjj values. The exclusive signal is therefore obtained 
by a fit of the Rjj distribution to expected signal and 
background shapes generated by Monte Carlo (MC) sim- 
ulations. The background shape used is checked with an 
event sample of heavy quark flavor dijets, for which ex- 
clusive production is expected to be suppressed in LO 
QCD by the J z = selection rule of the hard scattered 
di-gluon system, where J z is the projection of the total 
angular momentum of the system along the beam direc- 
tion 0. 



III. DETECTOR 

The CDF II detector, shown schematically in Fig. O 
is described in detail elsewhere The detector com- 
ponents most relevant for this analysis are the charged 
particle tracking system, the central and plug calorime- 
ters, and a set of detectors instrumented in the forward 
pseudorapidity region. The CDF tracking system con- 
sists of a silicon vertex detector (SVX II) [l2j, composed 
of double-sided microstrip silicon sensors arranged in five 
cylindrical shells of radii between 2.5 and 10.6 cm, and 
an open-cell drift chamber (l3j of 96 layers organized 
in 8 supcrlayers with alternating structures of axial and 
±2° stereo readout within a radial range between 40 and 
137 cm. Surrounding the tracking detectors is a super- 
conducting solenoid, which provides a magnetic field of 
1.4 T. Calorimeters located outside the solenoid arephys- 
ically divided into a central calorimeter (CCAL) [14], [l5| , 
covering the pseudorapidity range < 1.1, and a plug 
calorimeter (PCAL) [l6j, covering the region 1.1 < \r/\ < 
3.6. These calorimeters are segmented into projective 
towers of granularity Arj x A<fi « 0.1 x 15°. 

The forward detectors [l7|, which extend the cover- 
age into the r\ region beyond 3.6, consist of the MiniPlug 
calorimeters (MPCAL) 18], the Beam Shower Counters 
(BSC), a Roman Pot Spectrometer (RPS), and a system 
of Cerenkov Luminosity Counters (CLC). The MiniPlug 
calorimeters, shown schematically in Fig. 01 are designed 
to measure the energy and lateral position of particles 
in the region 3.6 < |ry| < 5.2. They consist of alter- 
nating lead plates and liquid scintillator layers perpen- 
dicular to the beam, which are read out by wavelength 
shifting fibers that pass through holes drilled through 
the plates parallel to the beam direction. Each Mini- 
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FIG. 3: Schematic drawing (not to scale) of the CDF II detector. 
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FIG. 4: Schematic cross sectional view of one of the two for- 
ward MiniPlug Calorimeters installed in CDF II. 



Plug is 32 (1.3) radiation (interaction) lengths deep. The 
BSC are scintillation counters surrounding the beam pipe 
at three (four) different locations on the outgoing pro- 
ton (antiproton) side of the CDF II detector. Covering 
the range 5.4 < \r]\ < 5.9 is the BSC1 system, which 
is closest to the interaction point (IP) and is used for 
measuring beam losses and for triggering on events with 
forward rapidity gaps. Lead plates of thickness 1.7 ra- 
diation lengths precede each BSC1 counter to convert 7 
rays to e + e~ pairs to be detected by the scintillators. 
The RPS, located at ~ 57 m downstream in the antipro- 
ton beam direction, consists of three Roman pot stations, 
each containing a scintillation counter used for triggering 
on the p, and a scintillation fiber tracking detector for 
measuring the position and angle of the detected p. The 
CLC [19(, covering the range 3.7 < \rj\ < 4.7, which sub- 
stantially overlaps the MiniPlug coverage, are normally 
used in CDF to measure the number of inelastic pp colli- 
sions per bunch crossing and thereby the luminosity. In 



this analysis, they are also used to refine the rapidity gap 
definition by detecting charged particles that might pen- 
etrate a MiniPlug without interacting and thus produce 
too small a pulse height to be detected over the MiniPlug 
tower thresholds used. 

IV. DATA SAMPLES AND EVENT SELECTION 

Three data samples are used in this analysis, referred 
to as the DPE, SD, and non-diffractive (ND) event sam- 
ples. The exclusive signal is derived from the DPE event 
sample, while the SD and ND samples are used for eval- 
uating backgrounds. The total integrated luminosity of 
the DPE sample is 312.5 ± 18.7 pb" 1 . 

The following trigger definitions are used: 

J5: a single CCAL or PCAL calorimeter trigger tower 
of E T > 5 GeV. 

RPS: a triple coincidence among the three RPS trigger 
counters in time with a p gate. 

BSCl p : a BSC1 veto on the outgoing proton side. 
The three event samples were collected with the following 
triggers: 

ND=J5, SD=J5 • RPS, DPE=J5 • RPS •BSCLj. 

The DPE events, from which cross sections are calcu- 
lated, were sampled at a rate of one out of five events to 
accommodate the trigger bandwidth. In the above sam- 
ple definition, ND events include SD and DPE contribu- 
tions, and SD events include DPE ones. This results in 
a "contamination" of background distributions by signal 
events, which is taken into account in the data analysis. 

The selection cuts used in the data analysis include: 

VTX cut (ND, SD, and DPE): no more than one re- 
constructed primary vertex within \z\ < 60 cm, imposed 
to reduce the number of overlap events occurring during 
the same beam-beam crossing at the IP. 

RPST cut (SD and DPE): RPS trigger counter pulse 
height cut, imposed to reject "splash" triggers caused 
by particles hitting the beam pipe in the vicinity of the 
RPS and spraying the RPST counters with secondary 
particles. 

JET cut (ND, SD, and DPE): events are required 
to have at least two jets with transverse energy E^* > 
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10 GeV within \rj\ < 2.5. The transverse energy of a jet is 
defined as the sum E^ e = YliEi sin($j) of all calorimeter 
towers at polar angles 9i within the jet cone. Jets are 
reconstructed with the midpoint algorithm (20| , which is 
an improved iterative cone clustering algorithm, using a 
cone radius of 0.7 in rj-cj) space and based on calorimeter 
towers with E T above 100 MeV. The E T of a jet is defined 
as the sum of the Et values of the clustered calorimeter 
towers. The jet Et is corrected for the relative response 
of the calorimeters and for the absolute energy scale. 

The above selection cuts define the DPE data sample 
(DPE) and are summarized below: 

DPE sample: (5) 

J5 ■ RPS ■ BSClp ■ VTX ■ RPST ■ JET. 

The DPE data sample consists of 415 688 events. 

Backgrounds in the DPE event sample fall into two 
general categories: (i) SD dijet events, in which the 
BSClp requirement is fulfilled by a downward BSC\ p 
multiplicity fluctuation to zero, and (ii) overlaps between 
a ND J5 trigger and a RPS trigger provided by either a 
low mass soft SD event that has no reconstructed ver- 
tex or by a scattered beam halo or ND event particle. 
To reduce these backgrounds, two more requirements are 
imposed on the data: a large rapidity gap on the out- 
going proton direction, LRG P , and passing the £p cut 
defined below. 

LRG p : this requirement is implemented by de- 
manding zero multiplicities in MPCAL p and CLC p , 
Nmpcal = ^clc = 0, added to the trigger require- 
ment of BSClp = 0. The LRG P approximately covers the 
range of 3.6 < 77 < 5.9. This selection cut enriches the 
DPE event sample in exclusive events by removing non- 
exclusive backgrounds. Although there is a substantial 
overlap between the pseudorapidity regions covered by 
MPCALp and CLCp, the requirements of MPCALp = 
and CLCp = are nevertheless complementary, as the 
two systems detect hadrons and electromagnetic parti- 
cles with different efficiencies. 

cut: 0.01 < £p < 0.12. In the high instantaneous 
luminosity environment of Run II, multiple pp interac- 
tions occurring in the same beam-beam bunch crossing 
may result in overlap events consisting of a ND dijet event 
overlapped by a soft SD event with a leading p trigger- 
ing by the RPS. These events, which are a background 
to both DPE and SD dijet events, can be well separated 
from diffractively produced dijet events using the variable 
£p , defined as 

1 N tower 

where the sum is carried out over all calorimeter towers 
with E T > 100 MeV for CCAL and PCAL, and E T > 20 
MeV for MPCAL. The tower Et and 77 are measured with 
respect to the primary vertex position. The variable £p 
represents the fractional longitudinal momentum loss of 




FIG. 5: (a) distribution of DPE events passing the LRG P 
requirement (solid histogram), with the shaded area repre- 
senting events in the region 0.01 < &f < 0.12; the dashed 
histogram shows the £p distribution for ND events passing 
the same LRG P requirement and normalized to the solid his- 
togram in the plateau region of 0.22 < < 0.50; (b) MP- 
CAL hit multiplicity, N^ /p , vs. BSClp hit counter multiplic- 
ity, SC1 , in SD events with 0.01 < < 0.12. 



the p measured using calorimeter information. For events 
with a gap on the p side, £jf is calibrated by comparing 
data with Monte Carlo generated events. Calibrated 
values were found to be in good agreement with values 
of £p measured by the RPS, ^ PS - Cm t ne proton side 
where there is no RPS, is obtained from calorimeter 
information using Eq. [6] in which —77 in the exponent is 
changed to +rj and is calibrated using the MC technique 
that was validated by the comparison with RPS data on 
the p side. Fig. [5] (a) shows t^f distributions for events of 
the DPE event sample selected with the LRGp require- 
ment. The events in the peak at £jf ~ 0.05 are dominated 
by DPE dijets, while the broad peak around ~ 0.3 are 
residual overlap ND dijet events for which the LRGp is 
caused by downward multiplicity fluctuations. 

In this analysis, we use the DPE dominated events in 
the range 0.01 < < 0.12. The same £p requirement is 
used in selecting the SD event sample. Figure[5](b) shows 
the MPCAL hit multiplicity, N p Mpi vs. BSC1 hit counter 
multiplicity, Ng SC1 , for the SD event sample. The ma- 
jority of the events have 5 < Nf ip < 10 and Ng SC1 > 3, 
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FIG. 6: (a) Mean Et, (b) mean rj of the two leading jets, and (c) azimuthal angle difference between the two leading jets of 
Et > 10 GeV in IDPE (circles), SD (solid histograms), and ND (dashed histograms) dijet events. 



but there are also some events with Ng Scl — N^. IP = 0, 
which are due to DPE events in the SD event sample 
efficiently passing the BSC1 P trigger requirement. 

The above trigger and offline selection requirements 
define the inclusive DPE event sample (IDPE). 

IDPE sample: (7) 

DPE ■ LRG P - ^ . 

The IDPE sample contains 20 414 events. 

In Fig. [51 we compare distributions of the mean dijet 
transverse energy, E^, = {E 3 T etl + Ei T et2 )/2, mean pseu- 
dorapidity, rj* = {rjjeti +%'et2)/2, and azimuthal angle 
difference, A0 = \(fij e ti — 4>jet2\, for IDPE (points), SD 
(solid histogram), and ND (dashed histogram) events. 
All distributions are normalized to unit area. The IDPE, 
SD and ND distributions exhibit the following features: 
(a) the E?p distributions for IDPE, SD and ND events 
are similar at low Ej* , but reach larger Ei^ values for SD 
and ND events due to the higher c.m.s. energies of IP -p 
and pp collisions relative to IP- IP collisions; (b) the ND 
rj* distribution is symmetric about rj* = 0, as expected, 
and the DPE distribution is approximately symmetric as 
the jets are produced in collisions between two Pomerons 
of approximately equal momentum (due to the approx- 
imately equal gap size on the p and p sides), while the 
SD distribution is boosted toward positive rj* (outgoing 
p direction) due to the jets being produced in collisions 
between a proton carrying the beam momentum, p^ 1 and 
a Pomeron of much smaller momentum, £ppa; and (c) 
the jets are more back-to-back in SD than in ND events, 
and even more so in IDPE events, due to less gluon radi- 
ation being emitted in events where colorless Pomerons 
are exchanged. 

V. INCLUSIVE DPE DIJET PRODUCTION 

The cross section for inclusive DPE dijet production is 
obtained from the IDPE event sample using the expres- 
sion 

incl _ N DPE( l ~ F BC) ,„s 
& DPE - JJ~ t ' ^ > 



where Nj^ pE is the number of DPE dijet events corrected 
for losses due to multiple interactions and for smearing 
effects on E^ 1 due to the detector resolution, Fbg is 
the non-DPE background fraction, L is the integrated 
luminosity, and e is the total event selection efficiency in- 
cluding detector acceptance. Details are provided below. 



A. Non-DPE background events 

There are two sources of non-DPE background events 
in the IDPE event sample underneath the DPE peak at 
0.01 < < 0.12 shown in Fig. M (a): one due to ND 
dijet events and the other due to SD ones. In both cases, 
the LRGp requirement of N PSC1 = N^ p = Np ic = is 
satisfied by downward multiplicity fluctuations. 

Non-diffractive background. The ND back- 
ground is caused by the RPS being triggered either 
by a real antiproton from an overlapping soft SD 
event or by a particle originating in beam-pipe or 
beam-gas interactions. This background is estimated 
from the distribution of ND dijet events with 

^bsci = Njfp = ^clc = normalized to the dis- 
tribution of DPE events in the region 0.22 < £f < 0.50, 
which is dominated by ND events. The DPE (normalized 
ND) £p distribution is shown in Fig. [5] (a) as a solid 
(dashed) histogram. Integrating the ND distribution 
over the range 0.01 < & < 0.12, we obtain the fraction 
of ND dijet background in the IDPE event sample to be 
Fgg = 13.3 ± 0.2 %. 

Single diffractive background. The SD back- 
ground is estimated by examining the correlation be- 
tween N psci + N^ /p and ^clc m * ne data sample. 
Figure[7](a) shows the distribution of N PSC1 +N^ /p ver- 
sus N P CLC for SD dijet events with 0.01 < £g < 0.12. 
The multiplicity along the diagonal, Ndiag, defined by 
^bsci + ^Mp=Np LC , is well fitted with a linear func- 
tion in the region 2 < Ndi ag < 14, as shown in Fig. 0(b). 
The diagonal distribution is used because it monotoni- 
cally decreases as Ndiag — ► providing the least back- 
ground under the peak. Extrapolating the fit to the bin 
with N^ SC1 = N^ /p = ^clc = y ields a SD back- 
ground fraction of F 00 — 24 ± 4 %. After correcting 
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FIG. 7: (a) Sum of the BSC1 hit counter multiplicity and 
MPCAL hit multiplicity, N^ gci + N^ p , versus CLC hit mul- 
tiplicity, N^ iC , in SD events with dijets of E^ 1 ' 2 > 10 
GeV and 0.01 < (,p < 0.12; (b) multiplicity distribution 
along the diagonal line in the left plot, Ndiag, defined by 
Nssci + ^^mp = ^clci w ith the solid line representing a lin- 
ear fit in the region 2 < ~Ndi ag < 9, and the dashed line the 
extrapolation of the fit to the Ndiag = bin. 



for a ND content of 42% in the SD data, estimated by 
applying the method used in evaluating F^q , we ob- 
tain a single diffractive background fraction of F^q = 
F 00 x (1 - 0.42) = 14 ±3 %. 



B. Corrections for multiple interactions 

Multiple interactions in the same beam-beam crossing 
may produce additional events which overlap the DPE 
event and cause it to fail the event selection requirements 
by contributing extra event vertices and/or by spoiling 
the rapidity gap on the proton side. Corrections for DPE 
event losses due to multiple interactions are considered 
separately for overlapping events with one or more recon- 
structed vertices, and for overlapping events which do not 
have a vertex but nevertheless spoil the LRG p . The lat- 
ter also account for LRG P losses due to beam background 
and/or detector noise. 

Overlap events with a reconstructed vertex. 
The average number of inelastic pp interactions per bunch 
crossing is given by n, = Li ■ <Jinei/fo, where d is the 
instantaneous luminosity, Ui ne i the inelastic interaction 
cross section, and / the Tevatron bunch crossing fre- 
quency of 1.674 MHz. The average number of pp interac- 
tions which have a vertex, ftf tx , is obtained by replacing 
Omei W1 th , the cross section of pp interactions with 
a vertex. From Poisson statistics, the probability that 
no pp interaction producing a vertex occurs in a beam- 
beam bunch crossing is given by P(0) = cxp(— h v i tx ). 
The number of observed DPE events, Nope, corrected 
for losses due to multiple interactions that yield overlap 
events with a vertex, W£ r p E , is obtained by weighting 
every event by P(0) _1 and summing up over all DPE 
events: NgRb - £,=7" exp(nf *). 

The value of of^, which is needed for evaluating n\ tx , 
is obtained from an analysis of the fraction of events with 
one or more reconstructed vertices contained in a sam- 
ple of zero-bias events collected by triggering on beam- 



beam crossings during the same time period in which 
the DPE sample was taken. The zero-bias sample is 
split into small sub-samples corresponding to different 
time slots of data taking to account for changes in beam 
and detector conditions, and the fraction of events with 
> 1 vertex function of instantaneous luminosity 

for each sub-sample is fit to the expected fraction given 
by 1 - P(0) = 1 - exp(-A • vYfa/fo) with as a 

free parameter. The average value obtained from the 
fits is cr^gj = 30.3 ± 1.5 (syst) mb, where the uncer- 
tainty is evaluated from the variations observed among 
the different sub-samples. Using this value, we obtain 
W£ r P r E = 189 317 ± 1325 events for the IDPE sample. 
The ±1.5 mb uncertainty in cr^ei leads to an uncertainty 
of ~ 2 % on N^p^, which is negligibly small compared 
to other uncertainties discussed below. 

Overlap events with no reconstructed vertex. 
The rapidity gap of DPE events remaining after rejecting 
events with more than one vertex could be further spoiled 
by the presence of additional soft pp interactions with no 
reconstructed vertex, by beam background, or by detec- 
tor noise. The correction for these effects is obtained 
from the same zero-bias samples by selecting events with 
no reconstructed vertex and evaluating the fraction F gap 
of events with LRG p . The correction factor, F~ ap , eval- 
uated for bins of different instantaneous luminosity and 
data taking time, is then applied to N c £ r P r E for the same 
instantaneous luminosity and time bins. Within the in- 
stantaneous luminosity range of 10 31 < Ci < 4 x 10 31 
cm _2 s _1 of our DPE data sample, F gap varies between 
70 % and 30 %. 



C. Event selection efficiency 

Jet selection efficiency. The trigger efficiency for 
jets with a calorimeter trigger tower of Et > 5 GeV is 
obtained from a sample of minimum-bias (MB) events 
triggered only on a CLC coincidence between the two 
sides of the detector. The E T e and rf et are reconstructed 
using the same algorithm as that used in the analysis of 
the IDPE dijet event sample. For MB events that contain 
a calorimeter trigger tower of E^ wer > 5 GeV, jets are 
selected if the trigger tower is contained within the Ar]- 
Acp cone of the jets. The trigger efficiency per jet is 
determined in bins of and rf et as the ratio of the 
number of jets containing a trigger tower of E!j? wer > 5 
GeV to the total number of jets in all MB events. The 
single tower trigger efficiency for a given DPE dijet event, 
esT, is derived from the efficiency per jet, the number of 
jets in the event, and the Et and 77 values of each jet. 
The DPE data are corrected for the trigger efficiency by 
assigning a weight of e^ T to each event. 

RPS trigger efficiency. The efficiency of triggering 
on a leading antiproton in the RPS trigger counters may 
be expressed as the product of the trigger counter ac- 
ceptance, Arps, and the RPS detector efficiency, epps- 
The latter can be further factorized into two terms: the 
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efficiency for finding the antiproton hits, eppsh, and the 
efficiency of the hit signals passing the trigger require- 
ment, ejipst- From a study of trigger counter signals 
produced by particles reconstructed as single tracks us- 
ing a zero-bias event sample, we obtain eppsh = 97±1 %. 
Using zero-bias events with signals in the trigger coun- 
ters consistent with the response expected for minimum 
ionizing particles, CRpst is found to be unity. The trig- 
ger counter acceptance is obtained from a simulation of 
SD events using the beam transport matrix to carry the 
recoil p from the IP to the RPS detectors. The total 
RPS acceptance for DPE dijet events is obtained from 
the expression 



a total 
A B.PS 
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where tp PS is the four momentum transfer squared mea- 
sured by the RPS and ^dpe the total number of DPE 
dijet events. For the events collected in our data taking 
period we obtain A^p^ = 78.4 ± 0.3 (stat.) %. 

£p cut efficiency. The requirement of 0.01 < £p < 
0.12 is used as a pre-selection cut to reduce ND dijet 
background due to superimposed pp interactions. How- 
ever, this cut also removes some DPE events. The ef- 
ficiency for DPE events retained by this requirement is 
obtained from the £g distributions of the DPE and ND 
dijet events shown in Fig. [5] (a) and used in Sec. IV Al 
to estimate the ND dijet background fraction in the 
IDPE data. Subtracting the normalized ND from the 
DPE events and evaluating the ratio of events within 
0.01 < £jf < 0.12 to the total number of events yields an 
efficiency of 98.5 ± 0.2 %. The deficit of this efficiency 
with respect to unity is due to fluctuations and calorime- 
ter resolution effects causing a small fraction of events to 
migrate outside the selected £^ range. 

Single vertex cut efficiency. The single vertex re- 
quirement (VTX cut), which is imposed to reject events 
with multiple interactions, also rejects single interaction 
events with extra misidentified vertices resulting from 
ambiguities in track reconstruction. Comparing the num- 
ber of IDPE events before and after imposing this re- 
quirement, we obtain a single vertex cut efficiency of 
eivtx = 98 ± 1 %. Using a similar method, the efficiency 
of the requirement of the vertex position being within 
\z\ < 60 cm is determined to be e zvtx = 92 ± 2 %. 

Jet reconstruction efficiency. The results pre- 
sented are based on events with at least two jets of E T e > 
10 GeV. The reconstruction of such relatively low Et jets 
in the CDF II calorimeters is prone to inefficiencies asso- 
ciated with the calorimeter measurement of particle ener- 
gies and the jet reconstruction algorithm used. Jet recon- 
struction efficiencies are studied using Monte Carlo dijet 
event samples generated with pythia 6.216 [2l[ and pro- 
cessed through a GEANT-based detector simulation [22l |. 
Simulated jets are reconstructed at both particle and 
calorimeter levels using the same jet reconstruction al- 



gorithm as that used in the data analysis. Then, events 
with matched pairs of particle and calorimeter level jets 
in y-(j> space are selected satisfying the requirement of 

Ai? = [{ycAL ~ VHAD? + (cbcAL ~ ^HAD?] 1/2 < 0.7, 

where y C AL (vhad) and (f> C AL {^had) are the rapidity 
and azimuthal angle of a calorimeter (particle) level jet. 
If more than one calorimeter level jet matches a hadron 
level jet, the closest matched calorimeter level jet is cho- 
sen. Using this method, the jet reconstruction efficiency 
ej e t, defined as the fraction of hadron level jets that have 
a matched calorimeter level jet, is obtained as a function 
of hadron level jet Et and ?y. We find that the value 
of 83 % at E 3 ^ = 10 GeV and reaches full 

jet 



efficiency at E T et ~ 25 GeV. The dijet reconstruction ef- 



ficiency for a given DPE event, edijet, is determined from 
the jet reconstruction efficiencies for the Et and r\ of the 
jets in the event. In evaluating cross sections, each DPE 
dijet event is assigned a weight of e^,- et , and the number 
of DPE dijet events is recalculated. 



D. Jet Et energy smearing 

The reconstruction of low Et jets suffers from energy 
smearing effects due to large fluctuations in the calorime- 
ter response to low Et particles. These effects, convo- 
luted with a steeply falling E T et spectrum, cause migra- 



tion of jets into adjacent E 3 T et bins. The smearing is 
unfolded as a function of E 3 T et using correction factors 
derived from inclusive DPE dijet even ts g enerated with 
the POMWIG Monte Carlo simulation [23|], described in 
Sec. IVI Al followed by a simulation of the detector. The 
E T et spectra of the second highest Et jet at particle and 
calorimeter levels are then compared. No matching be- 
tween particle and calorimeter level jets in y-<f> space is 
performed. The second highest E 3 T et is used in order to 
conform with the minimum E J T et 



thresholds imposed on 
in our cross section measurements and in available 
theoretical predictions. Calorimeter level jets are cor- 
rected for the relative response of the calorimeters and 
for the absolute energy scale. The correction factors, ob- 
tained for each E T et2 bin as the ratio of the number of 
particle level jets to the number of calorimeter level jets, 
vary from 0.93 ± 0.03 to 1.03 ± 0.03 within the region of 
10 < E T et < 50 GeV. This correction is applied to the 
measured inclusive DPE dijet cross section as a function 
ofE T et2 . 



VI. EXCLUSIVE DIJET PRODUCTION 

The exclusive dijet signal contained in the IDPE data 
sample is evaluated from the distribution of the dijet 
mass fraction, Rjj (=Mjj/Mx), by measuring the excess 
of events at high Rjj over expectations from the POMWIG 
Monte Carlo DPE event generator (23|, which does not 
simulate the exclusive process. Below, in Sec. IVI Al we 
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TABLE I: Diffractive/Pomeron structure function (DSF) used in the inclusive dijet pomwig Monte Carlo simulations. 



DSF Definition 



CDF0H1 


FH{P, Q' A ) of Eq. (H0|) for P p{p) and HI LO QCD fit 2 for P p(p) [2, 24] 


CDF 


F£ (p, g 2 ) of Eq. (H0|) for both P p and P p [24] 


Hl-fit2 


HI LO QCD fit2 with extended Q 2 range [28] 


ZEUS-LPS 


NLO QCD fit to ZEUS LPS data [30] 



demonstrate that the IDPE dijet data are well described 
by a combination of an inclusive MC generated distri- 
bution plus a non-DPE background obtained from data, 
in Sec. IVIBI we present the search for an exclusive dijet 
signal at high Rjj, in Sec. IVI Cl we discuss expectations 
from an exclusive dijet Monte Carlo simulation, and in 
Sec. IVI Dl we compare the data with an appropriately 
normalized combination of inclusive plus exclusive MC 
generated events. Cross section results for exclusive dijet 
production are presented in Sec. IVIII 



A. Inclusive pomwig Monte Carlo Simulation 

We first compare data distribution shapes with 
pomwig predictions to verify that the data are well de- 
scribed by the MC simulation apart from deviations ex- 
pected from the possible presence in the data of an exclu- 
sive dijet signal. The data used are the IDPE event sam- 
ple defined in Eq. (7) in Sec. IIV) which contains 20 414 
events. While this sample should contain a larger fraction 
of exclusive dijet events than the total DPE event sample 
defined in Eq. (5), it is used because in searching for an 
exclusive signal, agreement between pomwig predictions 
and data is more relevant if checked in a kinematic region 
as close as possible to that where the exclusive signal is 
expected. 

Dijet events are generated in pomwig using a 2 — > 2 
processes with Pomeron remnants (see Eq. [3]) and a min- 
imum transverse momentum cut of pJp m = 7 GeV/c. 
Each event is processed through the detector simulation 
and is required to pass the data analysis cuts. In com- 
parisons with IDPE data, the SD and ND backgrounds 
expected in the data are normalized to their respective 
14.0 % and 13.3 % values, estimated as described above in 
Sec. El and are added to the pomwig generated events. 
The MC distributions of pomwig DPE plus SD and ND 
background events and the corresponding data distribu- 
tions are normalized to the same area. Background SD 
distribution shapes are obtained from data satisfying the 
IDPE event sample requirements except for BSC P , which 
is replaced by N I BSC1 + N^ f p < 1 and Np ic < 1 exclud- 
ing events with N BSC1 = N^ fP = N^ LC = 0; ND shapes 
are extracted from J5 data satisfying the LRG P and 
requirements. 

As a diffractive/Pomeron structure function we use 
F?(/3,Q 2 ) oc p- 1 [H, where (3 is the longitudinal mo- 



mentum fraction of the parton in the Pomeron related 
to the x-Bjorken variable x Bj (x- value of the struck par- 
ton) by P = Xsj/£- The Q 2 dependence of F? is im- 
plemented as a weight to Fjj(P), determined from the 
CTEQ6L j2f| proton parton distribution function (PDF) 
at the Q 2 scale of the event. The justification for using 
the proton PDF is based on a Run II CDF measurement 
of a rather flat Q 2 dependence of the ratio of SD to ND 
structure functions, indicating that the Pomeron evolves 
with Q 2 similarly to the proton [26[ . We assign 46 % and 
54 % of the Pomeron momentum to quarks (u, d, u, d) 
and gluons, respectively, as measured by CDF in Run I 
from diffractive W, dijet, and 6-quark production (2?]]. In 
view of the above considerations, the following structure 
function form is employed in the MC program, 



0.46 • 



1 



E 



q—u,d,u.< 

+0.r>4- ~" r 



u q (Q 2 ) 



Pi 



UJg(Q 2 

Pa 



n 



(10) 



where uj q ^ g ){Q 2 ) is a weight used to include the Q 2 de- 
pendence of the quark (gluon) PDF and a = 10~ 5 is an 
arbitrary parameter employed to avoid a divergence at 
/3 = 0. 

Diffractive/Pomeron structure functions (DSFs) are 
also provided in the pomwig MC program, obtained 
from QCD analyses of HI diffractive DIS data 0. Two 
of the HI DSFs used in pomwig are the HI LO QCD 
fit2 (Hl-fit2) with a Q 2 range extended to 10 5 GeV 2 to 
cover the CDF range [H|, and the HI NLO QCD fit3 
(Hl-fit3). Recently, QCD analyses of diffractive struc- 
ture functions have also been performed by the ZEUS 
collaboration using diffractive DIS data obtained with a 
Leading Proton Spectrometer (LPS) [3fj| . and also by the 
rapidity gap (or Mx) method [31]. We have implemented 
programs returning NLO QCD fits for ZEUS-LPS and 
ZEUS-Mx structure functions for use in pomwig (see 
Ref. [12] for Mx data). However, a more recent QCD 
analysis of diffractive DIS data performed by HI using 
larger data samples and incorporating data from differ- 
ent final states [33| yields DSFs favoring the Hl-fit2 DSF 
over the Hl-fit3 DSF and in good shape agreement with 
the ZEUS-LPS DSF, while disfavoring the ZEUS-Mx 
DSF. Therefore, for consistency among measured DSFs 
at HERA, we exclude the Hl-fit3 and ZEUS-Mx DSFs 
from this analysis. 
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Guided by our Run I DPE dijet analysis results, in 
which Fjj measured from the ratio of DPE to SD dijet 
events was found to agree in shape and normalization 
with Hl-fit2, while the of Eq. (10]) measured from 
diffractive dijet production is suppressed by a factor of 
~ 10 relative to that from Hl-fit2, we use F? of Eq. (JTUJ) 
for the Pomeron emitted by the p (p) and Hl-fit2 for 
that emitted by the p (p) [26| ■ This combination, which 
will be referred to as CDFffiHl, is used as the default 
DSF in the POMWIG event generation. The four diffrac- 
tive/Pomeron structure functions used in the analysis are 
listed in Table UJ 

In Fig. [8] we compare (a) the average dijet E^* and (b) 
the average rf et distributions, E^. and 77*, between IDPE 
data and POMWIG generated events using CDF©H1, 
CDF, Hl-fit2, and ZEUS-LPS diffractive/Pomeron struc- 
ture functions. While all E^ distributions have simi- 
lar shapes, the data 77* distribution is broader than all 
simulated ones. The larger width of the data 77* distri- 
bution is due to the presence in the data of exclusive 
signal events concentrated in the pseudorapidity region 
around r\ ~ — 1 (see Sec. lVIDj) . Figure [9] shows data and 
POMWIG distributions of the dijet invariant mass, Mjj, 
and mass of the central hadronic system, Mx, where 

m x = [(EtT"^) 2 - (£^r er ^) 2 ] 1/2 > where E i 

is the energy of a tower with Et > 100 MeV for CCAL 
or PCAL and E T > 20 MeV for MPCAL and fk is a 
unit vector pointing to the center of the tower. Good 
agreement is observed between data and MC generated 
distributions for all four diffractive/Pomeron structure 
functions. 



B. Search for exclusive dijets 

We search for exclusive dijet production by comparing 
data with POMWIG simulated dijet mass fraction distri- 
butions, Rjj = Mjj/Mx, looking for an excess of data 
over simulation at high Rjj . Data and four pomwig Rjj 
distributions obtained with CDF©H1, CDF, Hl-fit2, and 
ZEUS-LPS DSFs are shown in Fig. [TO] All distributions 
are normalized to unit area. 

An excess of data over simulated events at high Rjj is 
observed for all four DSFs used in the simulation. This 
excess is examined for consistency with the presence in 
the data of an exclusive dijet signal by applying selection 
cuts expected to enhance the appearance of the signal. 
The following successive cuts have been studied: 

(a) LRGp : this cut, which is the equivalent of LRG P , 
retains events with N^ SC1 = N^ P = N^ LC . = 0, enforc- 
ing a gap approximately covering the region 3.6 < 77 < 
5.9. 

(b) E^ t3 < 5 GeV: third jet veto. Applying a veto 
on events with three (or more) jets of E^ t3 > 5 GeV 
further enhances the exclusive dijet signal, resulting in 
a narrower exclusive signal peak in the Rjj distribution. 
This requirement tends to shift events toward high Rjj 
values by removing, for instance, exclusive dijet events 



which contain extra reconstructed jets originating from 
gluon radiation in parton showers. In evaluating exclu- 
sive cross sections, the loss of such events is accounted for 
by correcting the data for the exclusive signal acceptance 
obtained from exclusive dijet MC simulations. 

(c) r)i et -cut: r? etx and/or rf et2 < -0.5. This cut 
exploits correlations in the t^' 1 vs. if et2 distribution, 
which is more symmetric around rf etx — 77 jef 2 = for 
inclusive than for exclusive events. 

Figure [TT] shows jf etl vs. -rf et2 distributions for DPE 
data satisfying all the above selection cuts, POMWIG 
dijet events generated with the CDF0H1 DSF, and ex- 
clusive dijet events generated with two different exclu- 
sive dijet MC simulations, ExHuME and ExclDPE, 
which are described below in Sec. IVI CI The pomwig 
generated events, which do not contain an explicit exclu- 
sive contribution, and the data, which are dominated by 
non-exclusive events, are scattered symmetrically around 
rjjeti _ ^jet2 _ ag requirements of LRG P and 
LRGp accept the same range of momentum loss frac- 
tions £p and £ p . Since, however, the recoil proton is not 
detected, the RPST trigger introduces a bias in the exclu- 
sive production case, resulting in an asymmetric Pp-P v 
collision with £p > £ p , boosting the dijet system toward 
negative 77. We exploit this kinematic effect by splitting 
the data and the events generated by each MC simulation 
into two samples, A and B, defined in if etl -rf et2 space as 
shown in Fig. [TT] The data samples A, for which at least 
one of the two leading jets has ff et < —0.5, contains most 
of the exclusive signal events, while sample B, comprising 
all other events, has a much reduced exclusive contribu- 
tion. In Fig. [T2] we compare IDPE data distributions 
of Ej., 77*, Mjj, and Mx for the background- rich region 
B with the corresponding pomwig distributions obtained 
using the CDF0H1 DSF. Reasonable agreement between 
data and simulation is observed. 



C. Exclusive dijet Monte Carlo simulations 

In the current analysis, we use two Monte Carlo event 
programs for generating exclusive dijet events: Ex- 
HuME 1.3.1 H and dpemc 2.5 [3|. ExHuME is 
a LO matrix element event generator founded on the 
perturbative calculations presented in Ref. [||, while ex- 
clusive dijet production in DPEMC 2.5 (ExclDPE) is 
based on the DPE non-perturbative Regge theory in- 
spired model of Ref. [|. In ExHuME, we generate ex- 
clusive events using the MRST2002 next-to-leading order 
proton PDF [36| . and implement parton showering and 
hadronization using pythia. In ExclDPE, we generate 
exclusive dijet events with the default parameters, us- 
ing HERWIG 6.505 [38| to simulate parton showering and 
hadronization. 

Distribution shapes of £!£,, 77*, and Rjj for POMWIG 
DPE events and for events generated by ExHuME and 
ExclDPE are compared in Fig.rT3] All distributions are 
produced using quantities reconstructed at the hadron 
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FIG. 8: (a) Mean transverse energy E't and (b) mean pseudorapidity 77* of the two highest Et jets in IDPE data (points) and 
POMWIG MC events (thick solid histograms) composed of POMWIG DPE signal (thin solid histograms) and the sum of SD and 
ND background events (dashed histograms). The data and POMWlG+background distributions are normalized to unit area. 
The POMWIG generated distributions in the plots (a) and (b) correspond to the four different diffractive/Pomeron structure 
functions used: CDF0H1, CDF, Hl-fit2, and ZEUS-LPS. 
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FIG. 9: (a) Dijet mass Mjj and (b) central hadronic system mass, Mx in IDPE data (points) and POMWIG MC events (thick 
solid histograms) composed of POMWIG generated DPE dijet events (thin solid histograms) and the sum of SD and ND generated 
background events (dashed histograms). The data and POMWlG+background distributions are normalized to unit area. The 
POMWIG DPE distributions in each set of four plots correspond to the four different diffractive/Pomeron structure functions 
used: CDF0H1, CDF, Hl-fit2, and ZEUS-LPS. 
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FIG. 10: Dijet mass fraction, Rjj, in IDPE data (points) and 
POMWIG MC events (upper histograms), composed of POMWIG 
DPE signal and the sum of SD and ND background events 
(lower dashed histogram) normalized to the background frac- 
tion in the data. The upper four histograms correspond to the 
four different diffractive/Pomeron structure functions used 
in POMWIG: CDF0H1 (solid), CDF (dashed), Hl-fit2 (dot- 
ted) and ZEUS-LPS (dot-dashed histogram). These four his- 
tograms and the data distribution are normalized to unit area. 



We first fit the Rjj distribution of inclusive DPE di- 
jet events satisfying the additional cuts (a) and (b) of 
Sec. IVIB1 but not requiring the f? jet -cut. Results arc 
shown in Fig. [T3J In plots (a) and (b) the two highest 
Et jets in an event are required to have Ej? ' > 10 GeV 
and in plots (c) and (d) > 25 GeV. The solid histogram in 
each plot is obtained from a binned maximum likelihood 
fit of the data with a combination of (i) POMWIG DPE 
plus SD and ND background events (dashed histograms) 
and (ii) exclusive signal events (shaded histograms) gen- 
erated by ExHuME for plots (a) and (c) or ExclDPE 
for plots (b) and (d) satisfying the same cuts as the data. 
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FIG. 11: Second jet rj vs. leading jet r\ for events with two 
jets of E J T etlJet2 > 10 GeV satisfying all IDPE requirements 



plus the additional requirements of LRGp and E T < 5 GeV 
(third jet veto): (a) data, (b) POMWIG generated events, (c) 
exclusive ExclDPE generated events, and (d) exclusive Ex- 
HuME generated events. The solid lines represent the ri Jet - 



level for events selected with E J T et1 ' 2 > 10 GeV, \r? etx > 2 \ < 
2.5, 0.03 < < 0.08, and 3.6 < \rf a v\ < 5.9. The 
Ej, spectrum is harder (much harder) in ExHuME (Ex- 
clDPE) than in POMWIG, while the dijet system in both 
ExHuME and ExclDPE is boosted toward negative ij* 
owing to the selected £p range, as explained above in 
Sec. IVTBl In the Rjj distribut ions, the exclusive jets 
emerge around Rjj ~ 0.8, while POMWIG events populate 1 
the low Rjj region. Both ExHuME and ExclDPE MC 
distributions exhibit a long tail extending toward small 10 
Rjj values due to gluon radiation. For comparison with 
data, the MC generated events are processed through a 1 
detector simulation. 



D. Comparison of data with combinations of 
inclusive and exclusive simulated events 



cuts of rf and/or rj J < 



-0.5 used to divide the rf -rf 



space into the exclusive signal-enriched region A and the non- 
exclusive dominated region B. 
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FIG. 12: Comparison of IDPE data distributions of (a) E^, 
(b) Tj* , (c) Mjj, and (d) Mx for the background-rich region B 
with corresponding pomwig distributions obtained using the 
CDF©H1 DSF; the events plotted are those used in plots (a) 
and (b) of Fig. ITT] 
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FIG. 13: (a) Mean of the two leading jets, (b) mean rf et , -q* and (c) dijet mass fraction, Rjj for pomwig DPE dijet 
events generated using the CDF©H1 DSF diffractive/Pomeron structure function (solid histograms), and for exclusive dijet 
events generated with the ExclDPE (dashed histograms) and ExHuME MC simulations (dotted histograms). All distributions 
are normalized to unit area. 



In each fit, the normalizations of the inclusive POMWIG 
and of the exclusive MC events are introduced as free pa- 
rameters. The Rjj data distribution is well reproduced 
within statistical uncertainties with both ExHuME and 
ExclDPE based exclusive contributions, yielding exclu- 
sive fractions of F excl — 15.0 ± 1.2 (stat.)% and F exc i — 
15.8 ± 1.3 (stat.)%, respectively, for E 3 T et1 ' 2 > 10 GeV. 

As a control check, we add the requirement of the 
?7? e *-cut and obtain distributions for data, pomwig, Ex- 
HuME, and pomwig+ExHuME events separately for 
samples B and A, shown in Figs. [15] (a) and (b), respec- 
tively. The relative normalizations of the total DPE data 
and MC event samples are fixed to those obtained in the 
fits of Fig. [TJ] (a), and the Rjj distributions are scaled 
by the number of events that pass the ry^'-cuts that de- 
fine the data samples. As expected, no significant ex- 
clusive contribution is observed in the data sample B. 
In sample A, good agreement is observed between the 
data and the pomwig+ExHuME combination, indicat- 
ing that the observed excess at high Rjj, whose frac- 
tion in the data has increased from 15.0 ± 1.2 (stat)% 
to 20.8 ± 0.8 (stat)% by the sample A selection cuts, is 
consistent with an exclusive signal in both shape and rel- 
ative normalization. Similar agreement is observed using 
the ExclDPE simulation. 

The fit of the data with MC simulated events would be 
expected to improve if the normalizations were left free 
to be determined by the fit. Binned maximum likelihood 
fits to the data of sample A are shown in Fig. [TB] The 
fraction of exclusive dijet signal, F exc i, is found to be 
23.0 ± 1.9 (stat)% for ExclDPE and 22.1 ± 1.8 (stat)% 
for ExHuME. 

In Fig. [TTl we compare the rj* distribution of the data 
of sample A with a MC generated distribution using 
POMWIG with the CDF©H1 diffractive structure function 
and an admixture of an exclusive signal of (a) 23% Ex- 
clDPE or (b) 22% ExHuME generated events, where 
the normalization was fixed to that obtained from the fits 
in Fig. 1161 Considering that the normalization was not al- 
lowed to very in performing the fit, reasonable agreement 



TABLE II: Fraction of exclusive dijet events in DPE dijet 
data, extracted from likelihood fits to data selected with 
the requirements of E 3 T et1 ' 2 > 10 GeV, E 3 T et3 < 5 GeV, 
rf et1 ' 2 > -2.5, rf 8 * 1 ^ < -0.5, 0.03 < < 0.08 and 
3.6 < |77 9ap j < 5.9, using combinations of pomwig+ExHuME 
or pomwig+ExclDPE distribution shapes. Results are listed 
for four different DSFs used in pomwig. The rj jetim < —0.5 
cut requires that at least one of the two highest Et jets be 
within rij e t < —0.5. Uncertainties are statistical only. 



DSF 


ExclDPE 


ExHuME 


CDF0H1 


23.0 ± 1.9 % 


22.1 ± 1.8 % 


CDF 


22.6 ± 1.9 % 


21.7 ± 1.8 % 


Hl-fit2 


26.0 ±2.1 % 


24.7 ± 2.0 % 


ZEUS-LPS 


25.4 ±2.1 % 


24.3 ± 2.0 % 



between data and simulation is observed, confirming our 
previous assertion that the broader data than pomwig 
simulated distribution seen in Fig. [5] is due to the exclu- 
sive contribution in the region around rf ~ —1. 

To determine the sensitivity to the DSFs used in the 
simulations, we have extracted the fraction F exc i using 
eight different combinations of DSFs, made up from each 
of the four DSFs used in pomwig (CDF©H1, CDF, Hi- 
nt 2 and ZEUS-LPS) with the DSF used in ExHuME 
or in ExclDPE. The eight F exc i values obtained, listed 
in Table [TTl are mutually consistent within the quoted 
statistical uncertainties. 



VII. RESULTS 

In this section we present results for both inclusive 
DPE dijet cross sections, o ,t ^p E , and for exclusive pro- 
duction, <Jjj Cl '. The inclusive cross sections are evalu- 
ated from the IDPE dijet data sample defined by the 
selection cuts listed in Eq. ([7]). Although the exclusive 
events are expected to be concentrated in the region of 
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F excl = 15.0 ±1.2% 
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FIG. 14: Dijet mass fraction in IDPE data (points) and best fit (solid histograms) with a mixture of (i) POMWIG generated 
events composed of POMWIG DPE signal and SD plus ND background events (dashed histogram), and (ii) exclusive dijet MC 
events (shaded histogram). The data and the MC events are selected from the respective IDPE samples after applying the 
additional veto cuts of LRGp and E J T et3 < 5 GeV. Plots (a) and (c) [(b) and (d)] present fits using ExHuME [ExclDPE] 
generated exclusive dijet events, while a requirement of £^? t2 > 10 GeV [25 GeV] is applied to plots (a) and (b) [(c) and (d)] . 



0.03 < £p < 0.08, as determined from simulations and 
from a sub-sample of the data with recorded RPS track- 
ing information, the larger range of 0.01 < < 0.12 
is used to ensure that there are no inefficiencies caused 
by resolution and/or possible systematic effects associ- 
ated with the calorimeter based definition of £p [26[ . The 
results are corrected for backgrounds falling within this 
larger £jf region. 

Exclusive cross sections are obtained by scaling cr D !0 E 
by the fraction of IDPE data that pass veto cuts (a), (b) 
and (c) of Sec. I VI Bl and multiplying the result by F exc i ■ 
A-exci-* wnere F exc i is the exclusive fraction and A exc i the 
acceptance of the veto cut(s) for exclusive events. As the 
veto cuts include a LRGp requirement, we recalculate the 
correction for spoiled gaps due to multiple pp interactions 
with no vertices. We then scale o^pe by the ratio of the 
correction for p®p gaps to that for only a p-gap and apply 
it in evaluating The acceptance of the exclusive 

cuts, A exc i, is obtained from the fraction of ExHuME or 
ExclDPE generated events passing the same cuts. 



In the following sections, we summarize the methods 
we used to calculate systematic uncertainties and their 
contributions to the total uncertainty. 



A. Jet Et smearing 

The large difference in the slope of the E^ 1 distribu- 
tions between inclusive POMWIG and exclusive MC gener- 
ated events seen in Fig. [TBI results in different corrections 
for smearing. Corrections for ExHuME and Ex- 
clDPE generated events are derived using the method 
described in Sec. IV Dl In obtaining the final results for 
Cjj Cl using ExHuME, the cross sections extracted by the 
above procedure are multiplied by the ratio of the correc- 
tions obtained from the ExHuME event sample to the 
POMWIG based corrections to account for the difference 
in Ej? spectra. 
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FIG. 15: Dijet mass fraction for IDPE data (points) 
and for POMWIG generated events (dashed histogram) com- 
posed of POMWIG DPE plus SD and ND background 
events, and for ExHuME generated exclusive dijet events 
(shaded histograms). The solid histogram is the sum of 
POMWlG©ExHuME events. Plot (a) shows distributions for 
event sample B, and plot (b) for event sample A. The events 
plotted pass all other selection cuts. The MC events are nor- 
malized using the results of the fits shown in Fig. [14] (a), 
scaled according to the actual number of events that pass the 
rf e -cut requirement. 



FIG. 16: Dijet mass fraction for IDPE data (points) and best 
fit (solid histogram) to the data obtained from a combination 
of POMWIG events (dashed histogram) composed of POMWIG 
DPE plus SD and ND background events, and exclusive dijet 
MC events (shaded histogram) generated using (a) ExclDPE 
or (b) ExHuME. The data and the MC events are from sam- 
ple A and are required to pass all other selection cuts. 



1. Jet energy scale 



B. Systematic uncertainties 



The systematic uncertainty in the exclusive fraction re- 
ceives contributions from uncertainties in the jet energy 
scale, unclustered calorimeter energy determination, jet 
trigger efficiency, jet Et smearing, non-DPE background, 
RPS acceptance, luminosity determination, knowledge of 
the diffractive structure function, statistics of MC event 
samples, underlying event determination, and the mod- 
eling of the underlying event. 



The uncertainty in E J T et associated with the jet en- 
ergy scale (JES) is evaluated by varying the uncertainties 
on the relative and absolute energy scale corrections by 
±lcr in estimating the efficiency for triggering on a single 
calorimeter tower of Et > 5 GeV, while simultaneously 
monitoring the number of jets with E T et above the de- 
sired threshold. Due to the steeply falling E 3 T et spectrum, 
the change in trigger efficiency increases with decreasing 
E 3 T et from - 2 3 l % for 10 < E%* < 15 GeV to +{\ % for 

25 < E T et < 35 GeV, resulting in a variation of the num- 
ber of IDPE dijet events accepted of ±21 % {tf 7 %) for 
E 3 T e 42 > 10 GeV (E T e 42 > 25 GeV). This is the dominant 
uncertainty in both the inclusive and exclusive dijet cross 
section measurements 
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data fits. The full difference between the number of inclu- 
sive events (exclusive signal fraction) obtained from the 
varied samples and that obtained from the default sam- 
ple is assigned as a systematic uncertainty on this cor- 
rection. The uncertainties on the exclusive cross sections 
propagated from these differences are ±13 % (±21 %) for 
E 3 T et1 ' 2 > 10 (25) GeV. 



3. Jet trigger efficiency 



+ Tf t2 )/2 



Jet trigger efficiencies have been discussed in Sec. IV CI 
The full difference of the efficiencies obtained from 
minimum-bias data and inclusive RPS triggered data is 
taken as a systematic uncertainty and propagated to an 
uncertainty on the exclusive signal fraction. 



4- Jet Et smearing 

Corrections for inclusive DPE dijets are obtained from 
samples of POMWIG MC dijet events. Statistical uncer- 
tainties on the correction factors are taken as systematic 
uncertainties and propagated to uncertainties on <J % ^p l E . 
The full difference between the exclusive dijet cross sec- 
tions obtained using corrections derived from ExHuME 
and ExclDPE MC generated events is assigned as a sys- 
tematic uncertainty on the exclusive cross sections. This 
uncertainty is ±4 % (±6 %) for E J T et1 ' 2 > 10 (25) GeV. 



FIG. 17: Comparison of mean pseudorapidity distributions 
rf between the IDPE data of sample A and the mixture of 
pomwig and (a) ExclDPE or (b) ExHuME simulated events 
normalized by the fit presented in Fig. [16] 



2. Unclustered calorimeter energy 

Uncertainties on the unclustered calorimeter energy 
scale affect the £p measurement, which in turn leads to 
uncertainties not only on the number of observed events, 
but also potentially on Rjj distribution shapes. However, 
the CCAL and PCAL energy scale uncertainties tend to 
cancel out in the Rjj ratio. 

Changing the energy scale of CCAL and PCAL by 
±5 % in calculating leads to a variation of ±1 % in the 

inclusive DPE dijet cross sections for £^f* ' > 10 GeV. 
Varying the MPCAL energy scale by ±30 %, indepen- 
dently of CCAL and PCAL, results in a cross section 
change of 7_[ %. This is due to more (less) events falling 
into the range of 0.01 < < 0.12 when the MPCAL 
energy scale is lowered (raised). The number of data 
events at high Rjj is less affected by these changes, since 
the MPCAL is less active for such events. After adjusting 
the jet energy scale, the fraction of exclusive dijet signal 
in the DPE data is re-evaluated by repeating the MC to 



5. Non-DPE background 

The dominant non-DPE background uncertainty is as- 
sociated with the SD background fraction of F§q = 
14 ± 3 %, which contributes an uncertainty on the cross 
sections of ±0.03/(1 - 0.14) = ±3.5 %. The ±3 % uncer- 
tainty on F§q has a negligible effect on the Rjj shapes 
relative to other uncertainties. 



6. RPS acceptance 

The acceptance of the RPS trigger counters, Apps, 
could vary with beam conditions and changing counter 
efficiencies. During the data taking period, Apps varied 
by at most ±6 %. This value is assigned as a systematic 
uncertainty on Apps and propagated to both inclusive 
and exclusive cross sections. 



7. Luminosity 

The luminosity uncertainty, which is applied to all 
cross sections, is 5.9 %, with 4.4 % due to the accep- 
tance and operation of the luminosity monitor and 4.0 % 
due to the uncertainty in normalization using the total 
pp cross section [39j. 
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8. Diffractive structure function 

We have examined the effect of the choice of DSF on 
comparisons between data and pomwig generated event 
distributions. Wc find that the Hl-fit2 and ZEUS-LPS 
DSFs yield similar kinematic distribution shapes, which 
are in reasonable agreement with the data, while the Hl- 
fit3 and ZEUS-Mx ones produce significantly different 
distributions and are clearly disfavored. Guided by these 
results, we use pomwig events generated with the Hl-fit2 
DSF to perform the fits to the data and evaluate F exc i 
and take the full difference between the exclusive cross 
sections obtained and those extracted using the default 
DSF of CDF0H1 as a systematic uncertainty. 



9. Statistics of Monte Carlo event samples 

Statistical uncertainties in the MC generated events 
are taken into account in MC to data likelihood fits per- 
formed to extract F exc i, so that the uncertainty in F exc i 
is due to the uncertainties in both data and MC event 
samples. The MC associated uncertainty is derived from 
the uncertainty in F exc i by quadratically subtracting the 
data uncertainty, determined from the number of events 
in the extracted signal, and propagated to the MC con- 
tribution to the exclusive cross section uncertainty. 



10. Underlying event 

The observed excess of data over simulated events at 
high Rjj in Fig.[TB]could be due to an overestimate of the 
underlying event (UE) activity in the simulation. We in- 
vestigated this possibility by following the methodology 
previously developed by CDF in generic UE studies in 
pp collisions [37|]. The rj-cj) space is split into three re- 
gions with respect to the leading jet axis, the "forward" 
(|A0| < 60°), the "transverse" (60° < Acj) < 120° or 
240° < A(f> < 300°), and the "away" region (120° < 
Acf> < 240°), and the UE is evaluated in the transverse 
region, which is sensitive to the particles outside the jets. 

Figure [TBI shows the number of calorimeter towers out- 
side the jet cones per event vs. tower detector ?y for IDPE 
data and the simulation. Good agreement is observed, 
except for a discrepancy at the highest \r]\ region, where 
the tower transverse size is smaller than that of hadron 
showers produced by particles interacting in the calorime- 
ter. This results in several tower "hits" per particle at 
high \r/\, which is difficult to accurately simulate. 

A more relevant UE comparison between data and 
simulation is that of transverse calorimeter tower mul- 
tiplicity (Nt) and Et distributions between data and 
MC generated events, as shown in Fig. [TS] for differ- 
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FIG. 18: The number of calorimeter towers per event with 
Et > 100 MeV which are outside the jet cones versus tower 
detector n for IDPE data (points) and pomwig MC events 
(thick line) composed of pomwig DPE signal (thin line) and 
the sum of SD and ND background events (dashed line) nor- 
malized to the background fraction per event. 



quantify the effect of this discrepancy on the extracted 
exclusive signal fraction, we compare transverse tower 
Et distributions between data and MC in the region of 
0.4 < Rjj < 0.7, which is in the plateau of the distribu- 
tion shown in Fig. [TB] (to avoid edge effects). Then, we 
modify the UE in the simulation to minimize the x 2 /d.o.f. 
and re-evaluate the exclusive fraction. Using the default 
MC simulation yields a x 2 /d.o.f = 1.4. The UE is mod- 
ified by scaling the tower Et of transverse calorimeter 
towers by a factor F(E T ) — Et X (1.5 - iV T /40) for 
Nt < 20, yielding a x 2 /d.o.f. = 0.4, which is lower by 
1 unit. The Et scaling has the effect of fewer calorime- 
ter towers being rejected by the calorimeter threshold 
cuts, thereby resulting in a larger exclusive signal frac- 
tion, F^ a c [ ed = 27.0 ±2.2 (stat)%. Since the discrepancy 
between the default and scaled MC distributions is only 
one unit of x 2 /d.o.f., we retain the default value for the 
fraction and use the scaled result to assign a systematic 



uncertainty of ±(F ( 



scaled 
exel 



F 



default 
cl 



default 
i: cl 



)/2 = ±9% 



11. Exclusive dijet model 

The full difference between the exclusive cross section 
values obtained using the ExHuME and ExclDPE MC 
ent Rjj bins in the range 0.5 < Rjj < 1.0. Agree- Rjj shapes is assigned as a systematic uncertainty asso- 



ment between data and simulation is observed, except 
in the low multiplicity and low Et regions where the 
data points fall below the MC generated histograms. To 



ciated with the exclusive signal modeling (see Table iDIf . 
This difference is mainly due to different amounts of ra- 
diation emitted from the jets. 
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TABLE III: Measured inclusive DPE and exclusive dijet cross sections, and the ratio of the exclusive to inclusive DPE dijet cross 
sections in the kinematic range E J T et1 ' 2 > Ef in , \rf etl ' 2 \ < 2.5, 0.03 < £ p < 0.08 (integrated over all t p ), and 3.6 < rj gap < 5.9. 



(Tope ± stat ± syst 



a jj C ' ^ stat ^ s V st 



ryexcl 



10 GeV 
15 GeV 
20 GeV 
25 GeV 
35 GeV 



14.5 ± O.ltgj nb 
1.43±0.02tg g nb 

267 ± 6t\T pb 
76.0 ± 2.7±||; e pb 

14.6 ± 1.2±|;| pb 



1.10±0.04to54 nb 
112 ± 7±% pb 
15.7±2.0^ 5 6 5 pb 



+ 4 ' 11 pb 
1.37±0.49ti;oi pb 



4.84±0.yo_ 328 

□+1.08 



7.6 ± 0.3±^ % 

7.8 ± 0.5t?'2 % 

5.9 ± OStil % 
6.4 ± l.3tf 6 9 % 
9.3 ± 3.4±^ % 
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FIG. 19: Calorimeter tower multiplicity and tower Et distri- 
butions in the transverse region with respect to the leading 
jet axis (60° < A</> < 120° or 240° < A0 < 300°) for four Rjj 
bins in the region of 0.5 < Rjj < 0.9. The tower multiplicity 
and tower Et distributions are normalized to unit area and 
to the number of transverse towers per event, respectively. 



in quadrature. The exclusive cross sections are plotted in 
Fig. [201 (a), where they are compared with hadron-level 
predictions of ExHuME and ExclDPE Monte Carlo 
simulations. The ExHuME predictions are favored by 
the data in both normalization and shape. The exclu- 
sive signal for E r et1 ' 2 > 10 GeV is established at a sig- 
nificance level of 6.1 er, determined from the value of 



r>excL 



7.6 ± 0.3±^ % presented in Table [TTIl The 



value of 6.1 a is obtained as the ratio of the central value 
of R%%°j = 7.6 divided by an uncertainty composed of 
the statistical uncertainty of 0.3% and the "downward" 
systematic uncertainty of 1.2%, combined in quadrature 
and yielding (0.3 2 + 1.2 2 ) 1 / 2 = 1.24 %, as the systematic 
uncertainty of -1.2% comes from an upward fluctuation 
of the background. 

In Fig. [201 (b), we compare the data exclusive cross 
section for events with Rjj > 0.8 plotted vs. jet E™ m 
with the ExHuME prediction and with the analytical 
calculation of exclusive dijet cross sections from Ref. [4l| 
(KMR). The <J e jj Cl is recalculated using the observed ex- 
clusive signal in the region of Rjj > 0.8. The KMR cross 
sections, which are based on a LO parton level calcula- 
tion of the process gg — > gg, have an 0(3) systematic un- 
certainty. The kinematic cuts used in KMR are slightly 
different from those used in the present analysis, which 
could lead to an effect of ~ ±20 % on the predicted cross 
section values [42[ . The good agreement seen in Fig. [201 
between the measured cr|J c ' and the KMR predictions 
multiplied by a factor of 1/3 suggests that the data are 
consistent with the KMR predictions within the quoted 
uncertainties. An even better agreement is reached by 
rescaling the parton transverse momentum in the KMR 
calculation to the measured jet transverse energy [43| . 

The ratio of exclusive to inclusive DPE dijet cross sec- 
tions measured from the data as a function of jet E™ m 
is shown in Fig. [201 (c). 



C. Cross sections 

Measured cross sections for inclusive DPE and exclu- 
sive dijet production, and the ratio of exclusive to inclu- 
sive DPE dijet cross sections for different E^f* 1,2 thresh- 
olds, are presented in Table IIIII The listed systematic 
uncertainties consist of all those discussed above added 



VIII. HEAVY FLAVOR QUARK JETS 

One of the most characteristic features of exclusive di- 
jet production is that at high dijet mass fraction it is 
dominated by the parton level process gg — > gg, as con- 
tributions from gg — > qq are suppressed. Born level cross 
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FIG. 20: Exclusive dijet cross sections for events with two 
jets of _Ey et > 10 GeV plotted vs. the minimum E^f 1 of 
the two jets in the kinematic range denoted in the figures: 
(a) total exclusive cross sections compared with ExHuME 
and ExclDPE predictions; (b) exclusive cross sections for 
events with Rjj > 0.8 compared with ExHuME (solid curve) 
and with the LO analytical calculation from Ref. [(J (see also 
Ref. [43|) scaled down by a factor of three (dashed lines) - the 
shaded area represents uncertainties in the calculation due to 
hadronization effects; and (c) the ratio of total exclusive to 
inclusive DPE cross sections. 



sections for exclusive production of a color-singlet dijet 
system of mass M are given by [44[ 



dt 



-{99 -» 99) 
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■Gw-flS) = 1 



9 na 2 
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(11) 
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where Et is the transverse energy of the final state parton 
and m q is the quark mass. The suppression of gg — ► qq is 
due to the factor (m 2 /M 2 )(l — Am 2 /M 2 ), which vanishes 
as m„ /M 2 -> (J z = selection rule [3). Exclusive 
gg — > qq contributions are also strongly suppressed in 
NLO and NNLO QCD, and in certain higher orders [irj |. 

The predicted exclusive gg-dijet suppression offers the 
opportunity of searching for an exclusive signal in IDPE 
data by comparing the inclusive dijet Rjj shape with that 
of data containing identified qq dijets. The presence of 
an exclusive dijet signal in the IDPE event sample would 
be expected to appear as a suppression in the ratio of qq 
to inclusive events at high Rjj . This data driven method 
avoids the use of MC simulations and can be used to cor- 
roborate the MC-based extraction of the exclusive signal 
from the inclusive data sample. As many systematic ef- 
fects cancel in measuring the ratio, a relatively small qq 
event sample can provide valuable information. 

To ensure quark origin, we select jets from heavy flavor 
(HF) b- or c-quarks, identified from secondary vertices 
produced from the decay of intermediate B or D mesons 
using the SVX II detector. Both b- and c-quark jets are 
used, since the suppression mechanism holds for all quark 
flavors. 

Below, in Sec. IVIIIA1 we describe the HF data sample 
and event selection requirements, in Sec. IVIII Bl we eval- 
uate the HF selection efficiencies and backgrounds, and 
in Sec. IVIlFCl we present the HF jet fraction results. 



A. Data sample and event selection 

The data used in this analysis were collected at a full 
rate (no pre-scaling) with a trigger satisfying the same re- 
quirements as the DPE trigger, Jet5+RPS+BSCl p , plus 
an additional one designed to enhance the HF jet content. 
The latter required the presence of at least one track with 
transverse momentum pt > 2 GeV/c displaced from the 
IP by a distance d of 0.1 < d < 1.0 mm, where d is the 
distance of closest approach of the track to the IP [4^ . 
The total integrated luminosity of this data sample is 
200±12pb" 1 . 

Jets are reconstructed using a CDF Run I based iter- 
ative cone algorithm [48| with an r\-(\> cone of radius 0.4. 
The SECVTX tagging algorithm is used to search for a 
displaced secondary vertex due to a B or D meson de- 
cay within a jet cone. This algorithm seeks tracks with 
hits in the SVX II within the jet cone, and reconstructs 
the secondary vertex from those which are significantly 
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FIG. 21: Displaced secondary vertex mass distribution for 
jets tagged by the SECVTX algorithm in DPE data (points). 
The solid histogram shows the result of a three-component 
binned maximum likelihood fit to the data, using Monte Carlo 
templates of distribution shapes, consisting of b-jets (shaded 
histogram), c-jets (dashed histogram), and other jets (dotted 
histogram) obtained from pythia dijet events. 



displaced from the primary vertex. A jet is considered 
SECVTX tagged if it has a secondary vertex consisting 
of at least two (or three) such tracks with pr > 1 (0.5) 
GeV/c. Events are further required to pass the IDPE 
selection criteria listed in Eq. ([7]). This selected "pretag" 
event sample contains 34 187 jets with at least two tracks 
in the SVX II. Applying the SECVTX tagging algorithm 
to the jets in the pretag sample yields 1,118 tagged jets 
with E 3 T et > 10 GeV and \rf et \ < 1.5. 



the data sample to which the algorithm is applied, and is 
therefore evaluated using a combination of a Monte Carlo 
simulation and HF jet fractions obtained from the data. 
The tagging efficiency for a b (c)-jet, €b(e c ), is obtained 
using pythia MC dijet events passed through a detector 
simulation, and is corrected for discrepancies observed 
between MC events and data. The fraction of b (c)-jets 
in the DPE data sample, Fb(F c ), is obtained from the 
fit to the M svtx distribution shape (Fig. [2Tj) . The tag- 
ging efficiency for a HF jet, efjg = F HF /[F b /e b + F c /e c ], 

is found to be 7.9 ± 1.4 % for E 3 T et ~ 15 GeV, where 
Fhf = Fb + F c . 

The data are corrected for efficiencies associated with 
the displaced track (DT) requirement in the trigger. 
These efficiencies are obtained from a data sample col- 
lected without requiring a displaced track. Selecting from 
this sample events that pass the DT requirement, we ob- 
tain the DT trigger efficiency for inclusive jets as the ratio 
e aiT = N^if /N a u, where N®f is the number of jets in the 
selected events and N a u is the total number of jets in the 
sample. Similarly, the efficiency for HF jets is obtained 
as e D H T F = NE T F /N HF = [N£J(1 - F 7 D J stag )]/[N tag {l - 
Fmistag)], where N tag (Nfi£) is the number of tagged jets 
and F mlstag {F^J stag ) the mistag background fraction in 



events without (with) the DT requirement. F t 



DT 

mistag 



evaluated from a 3-component MC template fit to the 
M svtx data distribution, consisting of b-, o, and other- 
jets at experimentally measured proportions (Fig. l2"Tj) . 



while F, 



mistag 



is measured from the corresponding M s 



fit to the DPE data collected without the DT require- 
ment. 



C. Heavy flavor jet fraction results 



B. Heavy flavor selection efficiencies 

A SECVTX tag in a jet without a HF quark is labeled 
as a "mistag." The mistag probability per jet, measured 
from inclusive jet data, is parameterized as a function of 
the number of tracks, Et, r], and <f> of the jet, and the 
sum over the Et values of all jets with E 3 T et > 10 GeV 
and \rj jet \ < 2.4 in the event. The mistag background in 
tagged jets is estimated by weighting each jet in the pre- 
tag sample by the mistag probability and summing up 
the weights over all jets in the sample. The total num- 
ber of mistag jets, N mistaff , is measured to be 104 ± 15. 
This number is consistent with a background estimate 
obtained from studies of the distribution of the invariant 
mass distribution M svtx of charged particles associated 
with a displaced secondary vertex (Fig. l2"Tj) . The mistag 
background for a given Rjj interval is evaluated by ap- 
plying the mistag probability to the jets in that interval 
of the pretag sample. 

The efficiency for tagging HF jets by the SECVTX al- 
gorithm depends on the composition of b- and c-jets of 



Results for the fraction F HF / incl of HF jets to all inclu- 
sive jets of E 3 T et > 10 GeV and \rji et \ < 1.5 for the IDPE 
event sample are shown in Fig. [2U (a) as a function of 
dijet mass fraction Rjj. The fraction is normalized to 
the mean value of the ratio of the HF to inclusive events 
over the four Rjj bins in the region of Rjj < 0.4, so 
that systematic uncertainties correlated among Rjj bins 
cancel out, as e.g. the uncertainties from corrections for 
data to MC tagging efficiency discrepancies or from the 
mistag background fraction estimate before and after the 
DT trigger requirement. Thus, 

—incl 

F H F/incl = {F}JF /incl) I -Rjj <0.4 ' &mcl ~ , (13) 

where <j mcl is the inclusive DPE jet production cross 
section only, o~ excl is the exclusive cross section, and 
{FHF/incl)\R sj <QA is the mean value of F HF/inc i in the 
range Rjj < 0.4 (45|. An exclusive dijet production rate 
contributing to the total rate but which is suppressed 
in HF dijet production would be expected to appear as 
a suppression in the fraction FnF/inci a t high Rjj [44| . 
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FIG. 22: (a) Measured ratio F HF / inc i of heavy flavor jets to all inclusive jets of > 10 GeV and \rj J \ < 1.5 as a function 
of dijet mass fraction Rjj, normalized to the weighted average value in the region of Rjj < 0.4, with systematic uncertainties 
represented by the shaded band; (b) values of F exc i^HF = 1 — F\ (filled circles) and F exc l^MC = 1 — i"a (open squares) as 



a function of Rjj, where Fi = F H 



F/incl 



I < F 



HF/incl 



> 



i, plotted on left vs. Rjj, and F2 is the ratio of POMWIG MC 



to inclusive dijet events obtained from the studies presented in Sec. I VII 
(F exc i^ H F systematic) uncertainties. 



the error bars (shaded band) represent statistical 



The suppression seen in Fig. [22] (a) is examined here for 
consistency with this hypothesis. 

The statistical uncertainties shown in Fig. [22] are 
dominated by the low statistics HF event sample. Mea- 
suring the fraction of HF to inclusive dijet events has 
the advantage of reducing systematic uncertainties com- 
mon to both event samples. The fraction is corrected for 
mistag backgrounds, tagging efficiency for HF jets, and 
displaced track trigger efficiencies for inclusive and HF 
jets. 

The systematic uncertainties on the ratio 
Fh F/inci j shown in Fig. fJH as shaded areas, are due to 
the uncertainties associated with the background and ef- 
ficiency estimates. The mistag background uncertainty 
is evaluated from the uncertainties associated with the 
determination of the mistag probability and propagated 
to an uncertainty in the fraction. The HF-jet tagging ef- 
ficiency, derived from combinations of PYTHIA MC gen- 
erated events and data, has an uncertainty propagated 
from the statistical uncertainty of the MC generated 
event sample, uncertainties from the correction for dis- 
crepancies between the tagging efficiencies derived from 
MC and data, and uncertainties from the b- and c-jct 
fractions in the tagged jet data sample. The displaced 
track trigger efficiency has two sources of systematic un- 
certainty: the statistical uncertainty of the IDPE sample 
collected without the displaced track requirement, and 
the uncertainty on the mistag background fractions be- 
fore and after applying the mistag requirement to the 
sample. In addition to the above uncertainties, we assign 
a systematic uncertainty associated with an increasing 
trend of c- to 6-jet fraction found in pythia generated 
events, which could contribute to a decreasing HF-jet 



fraction with Rjj due to the tagging efficiency e c being 
lower than ejj. 

To examine the consistency between the MC based ex- 
tracted exclusive dijet fraction and the data based sup- 
pression of the exclusive HF to inclusive dijet production 
rates we compare in Fig. fJH (b) the Rjj residual distri- 
butions defined as F exc i^MC = 1 — [MCi nc i/Datai nc i], for 
which the excess is defined as the inclusive DPE events 
observed above the inclusive Rjj distribution (composed 
of POMWIG MC events with ND and SD backgrounds) 
normalized to the DPE data at Rjj < 0.4 (open squares), 
and F exc i^HF = 1 — [F H F/inci/(FHF/inci)|fl 3 j<o.4] (filled 
circles) . The absolute values and Rjj dependence of the 
F eX ci^HF points in the region of 0.4 < Rjj < 1.0 are 
consistent with those of F exc i^MC > supporting an inter- 
pretation of the observed F HF / inc i distribution as a man- 
ifestation of the suppression of HF quark jets in exclusive 
production. 



IX. EXCLUSIVE DUETS AND DIFFRACTIVE 
HIGGS PRODUCTION 

The search for Higgs bosons is one of the top priorities 
of the LHC experiments. While the main effort of both 
the ATLAS and CMS experimental plans is directed to- 
ward searches for inclusively produced Higgs bosons, an 
intense interest has developed toward exclusive Higgs bo- 
son production, p + p-^p + H + p 49] . The exclusive 
Higgs production channel presents several advantages, in- 
cluding: (i) it can provide events in an environment of 
suppressed QCD backgrounds for the main Higgs decay 
mode of H — > bj e t + bj et , due to the J z — selection 
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rule discussed in Sec. IVIII1 (ii) the Higgs mass can be 
measured accurately with the missing mass technique by 
detectin g an d measuring the momentum of the outgoing 
protons (51J, (iii) the spin-parity of the Higgs boson can 
be determined from the azimuthal angular correlations 
between the two outgoing protons, and (iv) the method 
is universally sensitive to all exclusive Higgs production 
mechanisms. 

Theoretical predictions for exclusive Higgs boson pro- 
duction cross sections range from ~ 200 fb [Bj] to 2-6 
fb 50] for a Higgs boson mass of ~ 120 GeV at \/s=14 
TeV at the LHC. However, since exclusive Higgs boson 
and exclusive dijet production proceed through similar 
diagrams, as illustrated in Figs. [T|b (with no Pomeron 
remnants) and Fig. [2] the models can be calibrated by 
comparing their predictions for exclusive dijet cross sec- 
tions with measured values at the Tevatron. Further- 
more, measured exclusive dijet cross sections at the Teva- 
tron may also be used to evaluate backgrounds to the 
process H — ► bb from exclusive gg dijet production with 
gluons misidentified as 6-quarks in 6-tagging, or from b- 
quarks produced by gluon splitting, g — > bb. 
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FIG. 23: ExHuME exclusive dijet differential cross section at 
the hadron level vs. dijet mass Mjj. The filled points show 
cross sections derived from the measured crJJ ci values shown in 
Fig. [2Sl (top) using the procedure described in the text. The 
vertical error bars on the points and the shaded band rep- 
resent statistical and systematic uncertainties, respectively, 
obtained by propagating the corresponding uncertainties to 



the measured values of crj 



The solid curve is the cross 



section predicted by ExHuME using the default settings. 



A. M 7 , distribution 



The measured exclusive dijet cross section presented 
in Fig. [20] vs. jet E™ ln is converted to a cross section vs. 
dijet mass Mjj using the ExHuME Monte Carlo simu- 
lation with Mjj reconstructed at the hadron level. From 
the measured values of crfj Cl for the £7|f* ' thresholds 
given in Table. IIIH we obtain the cross section for each 
of the following E 3 T et2 intervals; 10-15 GeV, 15-20 GeV, 
20-25 GeV, 25-35 GeV, and 35 GeV or higher. After 
applying a hadron level E 3 T et2 cut, the ExHuME Mjj 
distribution for each E^ t2 interval is normalized to the 
cross section for that interval. Summing up over all the 
normalized Mjj distributions yields the ExHuME-based 
exclusive dijet differential cross section as a function of 
Mjj, d(jjj Cl /dMjj. The values obtained are corrected 
for a possible bias caused by the minimum threshold re- 
quirement of E^ 12 > 10 GeV by comparing the Mjj 
distributions with and without the Ei^ 12 cut. The de- 
rived d(jjj Cl /dMjj distribution is shown for Mjj > 30 
GeV/c 2 in Fig. [23] (solid circles). This distribution 
falls slightly faster than the default ExHuME prediction 
(solid curve) , as one would expect from the fact that the 
measured <T|J cZ (£'™ m ) falls somewhat more steeply with 
EJp in than that of ExHuME (Fig[20]), but overall there 
is reasonable agreement. This result supports the Ex- 
HuME prediction, and thereby the perturbative QCD 
calculation of Ref. [6j on which ExHuME is based. 



B. Higgs boson cross section 

From the ExHuME resulting values of da^ cl /dMjj, 
we obtain a^ cl « 360 fb for the range 115 < Mjj < 
145 GeV/c 2 , which corresponds to a ±12 % mass window 
around Mjj = 130 GeV/c 2 for jets within the kinematic 
region defined by the cuts denoted in Fig. [23] For SM 
Higgs boson production at the Tevatron, perturbative 
calculations [5(| predict crff c/ ~ 0.2 fb with a factor of 2-3 
uncertainty for a Higgs boson mass of inn = 120 GeV/c 2 , 
which leads to a ratio of exclusive Higgs signal to dijet 
background of Rn/jj ~ 6 x 10~ 4 . This value is in agree- 
ment with the estimate of Ra/jj — 6 x 10 -4 given in 
Ref. for tuh = 120 GeV/c 2 using an experimental 
missing mass resolution of AM missing = 3 GeV/c 2 at 
the LHC, rendering support to the prediction of the SM 
Higgs exclusive production cross section of 3 fb (with a 
factor of 3 uncertainty) presented in Ref. [6j. Measure- 
ments of exclusive dijet production rates in the Higgs 
mass range at the LHC could further constrain ajf cl 
through Rjj/jj. 

Models of exclusive Higgs production may also be 
tested using measured cross sections for exclusive 77 pro- 
duction, p+p — > p + 77+p, a process similar to exclusive 
dijet production. In the model of Ref. the 77 pro- 
duction is represented by the diagrams of Fig. [5] in which 
"jet" is replaced by "7". A recent CDF measurement 
[531 ] yielded a cross section upper limit close to the pre- 
dicted value, providing further support for this exclusive 
production model. 
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X. SUMMARY AND CONCLUSION 

We have presented results from studies of dijet produc- 
tion in pp collisions at y/s = 1.96 TeV using events with 
a leading antiproton detected in a Roman Pot Spectrom- 
eter and a forward rapidity gap on the outgoing proton 
side, collected by the CDF II detector during Fermilab 
Tevatron Run II. These events, presumed to be produced 
by double Pomeron exchange (DPE) , were extracted from 
a data sample of integrated luminosity 310 pb" 1 . In par- 
ticular, we have demonstrated the presence of exclusively 
produced dijets, p + p — > p + dijet + p, by means of de- 
tailed studies of distributions of the dijet mass fraction 
Rjj, defined as the dijet mass divided by the DPE sys- 
tem mass. In comparisons of data Rjj distributions with 
inclusive pomwig [23j Monte Carlo simulations, we ob- 
serve an excess of events in the data over the Monte Carlo 
predictions at high Rjj, which is consistent in terms of 
kinematic distribution shapes with the presence of an ex- 
clusive dijet signal as modeled by the ExHuME [34| and 
exclusive DPE in dpemc [35[ Monte Carlo simulations. 
To facilitate comparison with theoretical predictions, the 
exclusive dijet cross section, cr|J ci , and the ratio of ex- 
clusive dijet to inclusive DPE dijet cross sections have 
been measured as a function of minimum Et threshold 
of the two leading jets in an event. The measured values 
of afj cl favor the ExHuME over the dpemc predictions, 
and are found to be consistent with predictions from per- 
turbative calculations presented in Ref. [(|. 

The Monte Carlo based extraction of the exclusive dijet 
signal is checked experimentally using a largely indepen- 
dent sample of heavy flavor 6-tagged jet events extracted 
from 200 pb _1 of DPE data collected with a special trig- 
ger requiring a track displaced from the interaction point. 
As exclusive dijet production from gg — > qq is predicted 
to be suppressed by the Jz — selection rule relative 
to production through gg — > gg, the ratio of identified 



heavy flavor quark jets to inclusive jets is expected to 
decrease at high Rjj. For jets of JSjf > 10 GeV, we 
observe a suppression of the ratio of heavy flavor jets to 
inclusive jets in the region of Rjj > 0.4, which is consis- 
tent in shape and magnitude with the expectation from 
the exclusive signal extracted by the MC based method 
based on Ref. p. 

The present results, representing the first observation 
of exclusive dijet production in high energy pp collisions, 
provide a benchmark template against which to calibrate 
theoretical calculations of exclusive Higgs boson produc- 
tion. The prospects for an observation of exclusive Higgs 
boson production at the LHC have been briefly discussed 
in light of our measured exclusive dijet cross sections. 

Acknowledgments 

We thank the Fermilab staff and the technical staffs 
of the participating institutions for their vital contribu- 
tions. This work was supported by the U.S. Department 
of Energy and National Science Foundation; the Italian 
Istituto Nazionale di Fisica Nucleare; the Ministry of 
Education, Culture, Sports, Science and Technology of 
Japan; the Natural Sciences and Engineering Research 
Council of Canada; the National Science Council of the 
Republic of China; the Swiss National Science Founda- 
tion; the A. P. Sloan Foundation; the Bundesministerium 
fur Bildung und Forschung, Germany; the Korean Sci- 
ence and Engineering Foundation and the Korean Re- 
search Foundation; the Science and Technology Facilities 
Council and the Royal Society, UK; the Institut National 
de Physique Nucleaire et Physique des Particules/CNRS; 
the Russian Foundation for Basic Research; the Comision 
Interministerial de Ciencia y Tecnologfa, Spain; the Eu- 
ropean Community's Human Potential Programme; the 
Slovak R&D Agency; and the Academy of Finland. 



[1] G. Ingelman and P. Schlein, Phys. Lett. B 152, 256 
(1985). 

[2] T. Affolder et al. (CDF Collaboration), Phys. Rev. Lett. 

85, 4215 (2000). 
[3] Rapidity, y = | In , and pseudorapidity, r\ = 

— In tan | , are used interchangeably for particles detected 
in the calorimeters, since in the kinematic range of inter- 
est in this analysis they are approximately equal. 
[4] While the exclusive DPE process is part of the inclusive, 
it has become traditional to refer to the processes with 
and without Pomeron remnants as "inclusive" and "ex- 
clusive," respectively. In this paper we will follow this 
tradition. 

[5] A. Bialas and P.V. Landshoff, Phys. Lett. B 256, 540 
(1991). 

[6] V. A. Khoze, A. D. Martin, and M. G. Ryskin, Eur. Phys. 

J. C 14, 525 (2000). 
[7] M. Boonekamp, R. Peschanski, and C. Royon, Phys. 



Lett. B 598, 243 (2004), and references therein. 
[8] We use a coordinate system with origin at the center of 
the CDF detector, z coordinate along the proton beam 
direction, and y coordinate pointing up; the x coordinate 
points away from the center of the accelerator ring. 
[9] A. Berera, Phys. Rev. D 62, 014015 (2000). 

[10] V. A. Khoze, A. D. Martin, and M. G. Ryskin, Eur. Phys. 
J. C 19, 477 (2001) [Erratum- ibid. C 20, 599 (2001)]. 

[11] D. Acosta et al. (CDF Collaboration), Phys. Rev. D 71, 
032001 (2005). 

[12] A. Sill et al, Nucl. Instrum. Methods A 447, 1 (2000). 
[13] A. Affolder et al, Nucl. Instrum. Methods A 526, 249 
(2004). 

[14] L. Balka et al, Nucl. Instrum. Methods A 267, 272 
(1988). 

[15] S. Bertolucci et al, Nucl. Instrum. Methods A 267, 301 
(1988). 

[16] M. Albrow et al, Nucl. Instrum. Methods A 480, 524 
(2002). 



27 



[17] K. Goulianos, "Diffraction at the Tevatron: CDF Re- 
sults", presented at DIFFRACTION 2006 - Interna- 
tional Workshop on Diffraction in High-Energy Physics, 
September 5-10 2006, Adamantas, Milos island, Greece, 
published in Proceedings of Science, PoS(DIFF2006)016. 
M. Gallinaro et al, Nucl. Instrum. Methods A 496, 333 
(2003). 

D. Acosta et al., Nucl. Instrum. Methods A 494, 57 
(2002). 

G.C Blazey et al, "Run II Jet Physics: Proceedings of 
the Run II QCD and Weak Boson Physics Workshop"; 



ArXi v:hep-ex/0 005012 
T. Sjostrand et al., Comput. Phys. Commun. 135, 238 
(2001). 

E. Gerchtein and M. Paulini, eConl C0303241, 
TUMT005 (2003). 

B. E. Cox and J. R. Forshaw, Comput. Phys. Commun. 
144, 104 (2002). 

T. Affolder et al. (CDF Collaboration), Phys. Rev. Lett. 
84, 5043 (2000). 

J. Pumplin, D.R. Stump, J.Huston, H.L. Lai, P. Nadol- 
sky, W.K. Tung, "New Generation of Parton Distribu- 
tions with Uncertainties from Global QCD Analysis, " J. 
High Energy Phys. 0207:012(2002). 

K. Goulianos, "Diffraction at the Tevatron: CDF Re- 
sults," Proceedings of Science PoS(DIFF2006)016. 

F. Abe et al (CDF Collaboration), Phys. Rev. Lett. 78, 
2698 (1997); 79, 2636 (1997); T. Affolder et al. (CDF 
Collaboration), Phys. Rev. Lett. 84, 000232 (2000). 

C. Adloff et al. (HI Collaboration), Z. Phys. C 76, 613 
(1997). 

The maximum value of Q 2 in the HI x2 fit3 Pomeron 
PDF is 75 GeV 2 and is used as a default setting in the 
POMWIG generator. 

S. Chekanov et al. (ZEUS Collaboration), Eur. Phys. J. 
C 38, 43 (2004). 

S. Chekanov et al. (ZEUS Collaboration), Nucl. Phys. B 
713, 3 (2005). 

M. Groys, A. Levy, and A. Proskuryakov, "Proceed- 
ings Part B of the HERA and the LHC: A Work- 
shop on the Implications of HERA and LHC Physics"; 
|arXiv:hep ex/0601013| 

A. Aktas e t al. (HI Collaboration) , Eur. Phys. J. C 48, 
715 (2006); |arXiv:hep -ex/0606004| 

J. Monk and A. Pilkington, Comput. Phys. Commun. 

175, 232 (2006); [arXiv:hep-ph/0502077| 

M. Boonekamp and T. Kucs, Comput. Phys. Commun. 

167, 217 (2005); |arXiv:hep-ph/0312273| 

A. D. Martin et al, Eur. Phys. J. C 28, 455 (2003). 

A. Affolder et al. (CDF Collaboration), Phys. Rev. 



D65:092002, 2002. 

[38] G. Corcella et al, 3. Hig h Energy Phys. 0101, 010 (2001); 
|arXiv:hep ph/0210213| 

[39] S. Klimenko, J. Konigsberg, and T. Liss, Report No. 
FERMILAB-FN-0741 (2003). 

[40] K. Goulianos, Phys. Rep. 101(3), 169-219 (1983). 

[41] V.A. Khoze, A.B. Kaidalov, A.D. Martin, M.G. Ryskin, 
W.J. Stirling, "Diffractive processes as a tool for search- 
ing for new physics", in Gribov Memorial Volume, Pro- 
ceedings of the Memorial Workshop Devoted to the 75th 
Birthday of V. N. Gribov, Budapest, Hungary 22 - 24 
May 2005, pp.129-144, published by Word Scientific, Yu 
L Dokshitzer (LPTHE Universites Paris- VI- VII, Paris, 
France) and P Lvai & J Nyri (MTA KFKI RMKI, Bu- 
dapest, Hungary) editors; arXiv:hep-ph/0507040vl. 

[42] V. A. Khoze and M. G. Ryskin, private communication. 

[43] V. A. Khoze, A. D. Martin, and M. G. Ryskin, "New 
Physics with Tagged Forward Protons at the LHC," 
larXiv:0705.2"314l [hep-ph] . 

[44] V. A. Khoze, A. D. Martin, and M. G. Ryskin, Eur. Phys. 
J. C 23, 311 (2002). 

[45] We assume that in the region of Rjj < 0.4 there is no HF 
contribution from exclusive jet production (see Fig. I14|l , 

j.x,„j. _ l incl incl I incl it? VI 

so that ohf/o- = a HF /a = {i> HF/incl)\R jj <0A , 

where a^ F l is the HF jet cross section obtained from 

inclusive DPE events. 
[46] V. A. Khoze, M. G. Ryskin, and W. J. Stirling, 
"On radiative QCD backgrounds to exclusive H — > 

bb production at the LHC and a photon collider", 

|arXiv:hep-ph/0607134| fl. 
[47] W. Ashmanskas et al, Nucl. Instrum. Methods, A 518, 

532 (2004). 

[48] F. Abe et al. (CDF Collaboration), Phys. Rev. D 45, 
1448 (1992). 

[49] FP420 R&D Project: |http : //www . f p420 . coiiij "FP420, 
An R&D proposal to investigate the feasibility of in- 
stalling proton tagging detectors in the J^20m region of 
the LHC", CERN-LHCC-2005-0025. 

[50] A. B. Kaidalov, V. A. Khoze, A. D. Martin, and 
M. G. Ryskin, Eur. Phys. J. C 31, 387 (2003); ibid. C 
33, 261 (2004). 

[51] M. G. Albrow and A. Rostovtsev, Searching for the 
Higgs at hadron colliders using the missing mass method, 
|ArXiv:hep-ph /0009336 

[52] V. A. Khoze, A. D. Martin, M. G. Ryskin, and W. J. Stir- 
ling, Eur. Phys. J. C 38, 475 (2005). 

[53] T. Aaltonen et al. (CDF Collaboration), Phys. Rev. Lett. 
99, 242002 (2007). 



